Electromagnon Resonances in the Optical Response of Spiral Magnets 
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The optical response of spiral magnets is studied, with special attention to its electromagnon 
features. We show that these features trace back to the resonant magnetoelectric response resulting 
from the spiral ordering (irrespective of any concomitant ferroelectricity) . This response, being 
magnetoelectric in nature, not always can be reduced to an effective electric permittivity. We argue 
that electromagnons in spiral magnets can produce, in addition to the observed peaks in the optical 
absorption of multiferroics, a (dynamically enhanced) optical rotation and a negative refractive 
index behavior. 



The strong interplay between magnetism and ferroelec- 
tricity observed in a new generation of ferroelectromag- 
nets (or multiferroics) has prompted a renewed interest 
on magnetoelectric (ME) phenomena. In TbMnOa, for 
example, the electric polarization can be flopped by ap- 
plying a magnetic field [H and, conversely, the chirality 
of its magnetic structure can be changed by applying an 
electric field Q . The dynamic counterpart of these cross- 
coupling effects is the existence of hybrid magnon-polar 
modes, i.e., the so-called electromagnons, which also have 
been observed in the form of absorption peaks in optical 
experiments Q . The relatively large magnitude of these 
ME effects makes this type of materials very attractive 
as novel memory elements, optical switches, etc. 

The ME response is known to be an important ingre- 
dient in the electrodynamics of conventional magneto- 
electrics 0,0,0, [3|- III Cr203, for example, this response 
alone suffices to produce a nonreciprocal optical rota- 
tion In this paper we provide a continuum medium 
description of the dynamical ME effect in spiral mag- 
nets. We show that, in contrast to the static case, the 
dynamic ME effect has to be described in terms of two 
ME response tensors. The reason lies in the different 
reaction to the external perturbations: the electric field 
is seen as a force, while the magnetic field as a torque. 
Moreover, the frequency dispersion of these ME tensors 
reveals both magnon and polar modes in the form of res- 
onances. These resonances are further shown to be re- 
sponsible for electromagnon features, as observed in the 
aforementioned optical experiments on multiferroics (see 
also d, [IjIl^jEl)- We also show that these features, be- 
ing genuinely ME in origin, not always can be reduced to 
an effective electric permittivity as in previous interpre- 
tations. We also discuss briefly the possibility of having a 
dynamically enhanced nonreciprocal optical rotation and 
a negative refractive index behavior due to such a genuine 
ME response. 

Let us begin by recalling that, at the static level, the 
most general linear response of a (homogeneous) medium 
to the electric and magnetic fields, E and H respectively. 



can be expressed by the constitutive relations [j, [5| : 

P = x'E + aH, (la) 
M = d'^E-hx™H. (lb) 

Here P and M represent the electric and magnetic po- 
larization, and X™ are the electric and magnetic sus- 
ceptibilities, and a is the ME tensor of the medium. The 
same tensor a enters in these two equations {ajj = aji) 
because it traces back to the same coupling —aijEiHj in 
the free energy of the system. Only a restricted number of 
magnetic symmetry classes allow for this linear coupling, 
in which case the system is termed as magnetoelectric. 

In the case of spiral magnets the inhomogeneous ME 
effect [13, EH [3] is always at work. This effect describes 
the (universal) coupling between the electric polarization 
and nonuniform distributions of magnetization. For our 
purposes it can be taken as [l5| 

-/P • [M(V -M) - (M • V)M]. (2) 

This coupling is believed to be behind the ferroelectric- 
ity of -RMnOa compounds, as comes from the observation 
that a cycloidal magnetization generates a term — P • E^s 
out of this coupling, where the constant Ecff reflects the 
lack of inversion symmetry in the cycloidal. This is the 
so-called spiral or spin-current mechanism for ferroelec- 
tricity, in which the polarization appears in the plane of 
the cycloidal perpendicular to its wavevector [l^. To 
describe more complicated cases one has to go beyond 
^ and consider different sublattice magnetization fields 
and/or the transformation properties of P and M under 
the elements of the corresponding point group (171]. 

To extend Eqs. ((T|) to the frequency domain we have 
to deal with the dynamics of the system. In this dynam- 
ics the coupling ([2]) already leads to the hybridization of 
magnons with polar modes (see below). It is likely that 
a more elaborated version of this coupling is needed to 
explain, e.g., the electromagnon selection rules observed 
in multiferroics for which the multisublattice magnetic 
structure seems to play a role Nevertheless we re- 
strict ourselves to the "minimal" coupling ^ because, 
already at this level, the dynamical ME effect can be 
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shown to be richer than noticed before. Compared to the 
static case, for instance, the dynamical ME effect exhibits 
the following asymmetry. Any hybridization implies that 
magnons can act as effective (time-dependent) electric 
fields for phonons and vice versa, phonons can act as ef- 
fective magnetic fields for magnons. Then charges will 
try to follow the mangon-induced field that is pushing 
them back and forth, whereas spins will tend to precess 
about the phonon-induced one (here we have a torque). 
This eventually translates into two dynamical ME ten- 
sors, in contrast to the static case ([T]). 

Let us compute these dynamical tensors. In the pres- 



ence of an electromagnetic radiation, electric and mag- 
netic polarizations will deviate with respect to the cor- 
responding background distributions: P = P^°^ -I- p and 
M = ]V['^°) -I- m. These deviations p and m are assumed 
to be small (proportional to the external fields), so the 
equations of motion for P and M can be linearized with 
respect to these quantities. For the sake of simplicity, 
we restrict ourselves to background magnetizations con- 
taining only one periodicity (i.e., with wavevectors ±Q). 
Thus, if we take the equation of motion for the electric 
polarization, in Fourier space we get 



J 



i(q,c.)p(q,c.)«E(q,c.) + 2z/ ^ [(q' • M(")(q'))m(q - q',c.) - M(0)(q')(q' • m(q - q',c.))] (3) 



q'=±Q 



r 



in the limit g <C Q. Here A represents the inverse electric 
susceptibility (in the absence of ME coupling A^^ = x"^)- 
In this equation we can see that, in fact, polar modes 
are linearly coupled with the deviations m(q ± Q,cj) of 
the magnetic structure by virtue of the modulation of this 
latter, i.e., we have electromagnons. Close to the magnon 
frequencies, as in ordinary antiferromagnets, the dynam- 
ics of the magnetization is expected to be described by 
the Landau-Lifshitz equation. Then, the nonlinear char- 
acter of this equation, together with the non-uniform 
magnetic background M*^°\ make it possible the linear 
coupling between these excitations and long-wavelength 



external fields (see e.g. [13|, |18[). Neglecting the ME 
coupling here for a while we have 

m(q±Q,c.) x'"-±Q(q,a;)H(q,^), (4) 

where the poles of ^"'^Q are associated with the char- 
acteristic excitations of the modulated structure 13, T§ |. 
Substituting this expression into ([3]) we obtain 



P(q, = X*'(q, t^)E(q, uj) + ci(q, w)H(q, lo), (5) 



where 
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a,j(q,cj)=2i/ q'kM[°\q'){SefSkk' -Sek'Skr)Xu'ifi,'^)xP} ''(q,^) 

q'=±Q 



(6) 



r 



The constitutive equation ^ replaces (jlal) for dynami- 
cal processes. The dynamical ME tensor a traces back 
to the magnon-phonon hybridization resulting from the 
ME coupling [Eq. ([2]) in our case]. Both these modes 
yield resonances in d, and these resonances are further 
responsible for electromagnon features in optical experi- 
ments (see below). 

As mentioned before, the fact that polarization and 
magnetization dynamics are different produces a certain 



asymmetry in the dynamical ME response. Carrying 
out similar manipulations the equation of motion for the 
magnetization can be reduced to an expression analogous 
to Hbl): 



(7) 



m(q, oj) = /3(q, w)E(q, uj) + x™(q, t^)H(q, uj), 
where 



J 



P^j{^_,uJ)=i-ffeu■^" Y qkAPj' iq,')M>"' i-c^'){^j'^,5kk' - (5j'fc'fe')Xi'j(q, ^) 



(8) 



q'=±Q 
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with 7 the gyromagnetic factor. As we see, this tensor 
(3 is not the mere transpose of the tensor a given in Q 
and, in contrast to d, does not contain information about 
magnon excitations (x"*'^*^ does not enter here). 

Let us now consider specific examples of magnetic 
structures. The first structure discovered with a long- 
period modulation was the helical one [lit : 

M^"^ (r) = Ml cos(Q • r) A + Ah sin(Q • r) z, (9) 

where Q = Q y. This type of magnetic ordering is ob- 
served, for example, in CaFeOs [20|. In this case, the 
inhomogeneous ME coupling ^ is ineffective in produc- 
ing an electric polarization since this structure does not 
break inversion symmetry. Nevertheless, it gives rise to 
a dynamical ME effect. To the lowest order (i.e., con- 
sidering the external field as the effective field acting on 
the magnetization in the Landau-Lifshitz equation), the 
non-zero components of the susceptibility are 



m,±Q ^ _ m,±Q 

m,±Q ^ _Y™.±Q rsr 



Mi°)(±Q), 

M(°)(±Q). 



(10a) 
(10b) 



If the electric susceptibility is diagonal, this means 
that the non-zero components of the dynamical ME ten- 
sors are a^x, ct^z, Pxx and f3zz- The ME response gener- 
ated dynamically in this case turns out to be analogous to 
the static one of the prototypical magnetoelectric Cr203 
(see, e.g., 0,0). 

Another important class of magnetic distributions is 
the cycloidal one: 

M(°)(r) = A/2Cos(Q-r)y + M3sin(Q-r)z, (11) 

with Q = Q y. The magnetization in RM11O3 com- 
pounds, for example, develops this type of modulation, 
and its appearance is accompanied with ferroelectricity 
as we have explained above. The dynamic ME response 
in this case has the following features. The susceptibility 
T^ni,±Q jtqj, cycloidal has the non-zero components: 



m,±Q ^ _ m,±Q 



A-xz 



Vyx 



Mi°)(±Q), 
Mf (±Q), 



(12a) 
(12b) 



to the lowest order. In consequence, for a diagonal elec- 
tric susceptibility, the only nonzero component of the ME 
tensor a turns out to be 

a,y{q,io) = 4z/QM2xL(q,^)x",''^(q,^)- (13) 

The fact that axy 7^ ayx{= 0) is a consequence of the 
inequivalence between the x and y directions in this mag- 
netic structure. This refiects also in /3, whose non-zero 
component reduces to (3yx as can be seen from (|8|). 

The optical response of /JMnOa compounds with this 
type of cycloidal magnetization shows absorption peaks 
at frequencies too small to be connected with pure 



phonon modes (- THz) [1, S 0, In TbMnOs, 

in particular, these peaks have been correlated with the 
low-lying excitations of the cycloidal observed by inelas- 
tic neutron scattering experiments [2l|. So they are in- 
terpreted as due to the electromagnon excitations natu- 
rally expected from the coupling ^ (see, e.g., [l!])- The 
electromagnon response to an external ac electric field 
is computed in [22{ as the (fluctuation) contribution to 
the electric permittivity due to the cycloidal excitations, 
and these results are further used to derive certain se- 
lection rules for the above optical experiments (see e.g. 
[lot). One has to realize that, however, this is not the 
whole story. Elcctromagnons actually react to both elec- 
tric and magnetic components of the external radiation 
due to their hybrid character, so their final response can 
be more complex [Eqs. ([5]) and ([7|), in general, do not 
reduce an effective permittivity] . 

Let us illustrate this point by computing the reflec- 
tion coefficient for a vacuum-cycloidal magnet interface. 
For the sake of concreteness we restrict ourselves to the 
case of normal incidence and linear polarization along the 
principal axes of the magnet [which are assumed to be 
the axes of the cycloidal (fTTj) in the following] . The result 
still depends on the orientation of the incident field with 
respect to the cycloidal. If the wavevector of the cycloidal 
Q is parallel to the interface the process is insensitive to 
the dynamical ME effect \a and (3 do not enter the re- 
flection coefficient, which is given by the standard Fresnel 
formula (see e.g. [5|])]. The same happens if Q is perpen- 
dicular to the interface and the electric field is along the 
polar axis of the cycloidal. However, if the electric field 
is perpendicular to the polar axis (i.e., the incident fields 
arc E* II X and H' || z), the ME effect comes into play. 
In this case, the dispersion law for the light propagating 
through the magnet is 



(14) 




and the reflection coefficient is found to be 



(15) 



(hereafter wc drop subindices since we are dealing with 
the only non-zero components of the ME tensors). In 
these expressions we see that the dynamic ME effect 
results in a new (effective) permittivity Scff = £ ~ ^ 
that now has poles at the magnon frequencies because of 
a ~ T^™'Q [see Eq. ([T^]. i.e., this effective permittivity 
has electromagnon features. This is in tune with (2^ and 
the general interpretation of the experimental data (see 
e-g- [HI)- For other orientations, however, the actual 
situation turns out to be a bit more subtle. 
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If. for example, the plane of the cycloidal is parallel 
to the interface and the electric field is directed along 
the polar axis (E* || z and H* || y), the ME couphng 
effectively result in the modification of the magnetic per- 
meability (not the electric permittivity as before). This 
is not captured in Q because only the influence of the 
electric field is taken into account. Experimentally, how- 
ever, no electromagnon feature seems to be observed for 
this orientation [lO| . The reason may be the fact that the 
magnetolelectric effect is here reduced just to an effective 
permeability (since the magnetic response alone is gener- 
ally quite weak) and/or the necessity of going beyond the 
isotropic coupling ^ for these materials. Furthermore, 
if the electric field is perpendicular to the plane of the 
cycloidal and this plane is perpendicular to the interface 
(E' II X and H* || y), the dispersion law is obtained from 
the equation: 



' ck 



(16) 



In this case the ME effect plays a genuine role, not re- 
ducible to a mere modification of the electric permittivity 
(or magnetic permeability) as before. The waves associ- 
ated with the two solutions of p6)) . for example, have 
different phase velocities. This possibility of removing 
the degeneracy between forward and backward waves is 
known long ago in genuine magnctoelectrics 0, 0] . 

It is worth mentioning that the field transmitted into 
the magnet can acquire a longitudinal component due 
to the ME effect. In the case E' || x and H' || z, for 
example, the transmitted field is such that 



Hi 



f3 



[e- 



af3\ 



(17) 



This possibility is also known for ordinary magncto- 
electrics. To probe experimentally this longitudinal com- 
ponent can be somewhat difhcult, but there is a related 
aspect of the ME effect whose experimental verification 
is, at least conceptually, much easier. It is the possibility 
of having a ME rotation of the reflected light. This possi- 
bility is quite obvious for the helical structure ^ taking 
into account that its ME response is analogous to that of 
Cr203 as we have seen. A similar rotation (see, e.g., 0) 
is therefore expected, with the particularity that in heli- 
cal magnets like CaFeOa it can be significantly enhanced 
due to the resonant behavior of the ME response. This 
dependence on the frequency is absent in conventional 
magnctoelectrics [l^]. In spiral magnets, it is closely re- 
lated to the electromagnon features in their spectrum 

Let us now explore yet another phenomenon that 
might benefit from these features. The ME effect has 
been pointed out as an interesting possibility to achieve 
a negative refractive index behavior [25l . [jq . recently 
demonstrated experimentally pTj . The key point in these 



experiments is the fabrication of metamaterials with chi- 
ral constituents. As in ordinary negative index metama- 
terials, the achievement of a negative index regime relies 
on the resonant response of the resulting system. In these 
cases, one basically deals with the resonances of the con- 
stituent particles [26|. This imposes severe limits to the 
range of frequencies at which the corresponding negative 
index behavior can be achieved. In the case of spiral mag- 
nets, on the contrary, it is the collective behavior of the 
system what gives rise to the ME effect. But the resulting 
resonant behavior [see Eqs. ^ and ([5|] is basically the 
same than in chiral metamaterials [26[. Consequently 
these type of magnetic structures may also result in a 
negative index behavior, now at the frequencies of the 
corresponding electromagnons (^^ THz for natural com- 
pounds). 

It is worth mentioning that spatial dispersion effects 
0, UHl can also be dynamically amplified in spiral mag- 
nets. Generally spatial dispersion produces minute ef- 
fects in optical experiments and, in practice, the response 
of the system is well described by the limiting g — > be- 
havior of the electric and magnetic susceptibilities 
0, [111 . Accordingly, in the computation of the ME re- 
sponse tensors a and $ we have neglected terms 0{q) 
coming from the inhomogeneous ME coupling ([2|) [sec Eq. 

Close to the electromagnon resonances, however, 
this neglection might be unjustified since these terms are 
dynamically enhanced by the same resonant mechanism 
that operates for ce{q = 0) and ${q = 0). The inhomoge- 
neous ME coupling then has to be considered in its full 
extent [l5| . and spatial dispersion effects may compete 
with the dynamic q = ME response in a similar way 
that it does, for example, with the static ME effect in 
CrsOa 0- 

In summary, we have shown that the optical response 
of spiral magnets exhibits electromagnon features en- 
coded in the form of a resonant magnetoelectric response. 
Spiral ordering does not have to be accompanied with 
multiferroicity (and/or a static magnetoelectric effect) to 
have these features. We have discussed the role of this dy- 
namical response in optical experiments on multiferroics, 
showing that the observed electromagnon features not al- 
ways can be reduced to an effective electric permittivity. 
We also have argued that electromagnon resonances in 
spiral magnets can amplify spatial dispersion and non- 
reciprocal effects. These resonances may also permit to 
achieve a negative refractive index behavior. 

I acknowledge P. Bruno, A. Levanyuk and especially 
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